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Abstract 

We study the existence of periodic solutions of the non-autonomous periodic Lyness ' 
recurrence u n+ 2 = (o n + Un+i)/u n , where {a n } n is a cycle with positive values a,b and 
with positive initial conditions. It is known that for a = b = 1 all the sequences 
generated by this recurrence are 5-periodic. We prove that for each pair (a, b) ^ (1, 1) 
there are infinitely many initial conditions giving rise to periodic sequences, and that 
the family of recurrences have almost all the even periods. If a ^ b, then any odd 
period, except 1, appears. 
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1 Introduction and main results 



The study of periodic orbits in periodic non-autonomous discrete dynamical systems is 
a classical topic that, mainly driven by some conjectures in mathematical biology, has 
attracted again the researcher's interest in the last years. This is because these ones are 
good models for describing the dynamics of biological systems under periodic fluctuations 
whether due to external disturbances or effects of seasonality (see [H [TUl EE21 E3 [23] for 
instance and references therein). On the other hand, the existence of discrete integrable 
systems with periodic coefficients is currently being the focus of some mathematical physics 
research both on those systems belonging to the celebrated QRT family of maps or 
on other families as the HKY one (see \17\ [22] and [16] , respectively) . Our work is in the 
interface of these two streams, since the main goal of the paper is to characterize the global 
periodic structure of a particular family of QRT maps: the one associated to the 2-periodic 
Lyness equations. 

Nowadays, the dynamics of the autonomous Lyness' difference equation 

a + u n +\ 

u-n+2 = with a > 0, and u±, U2 > 0. (1) 

u n 

is completely understood after the research done in [H [3] and [26] (see also |T5] ) . In summary, 
the dynamics of equation can be studied through the dynamics of the map F a (x,y) = 
(y, (a + y)/x). This map has a first integral V a such that their level sets in Q + := {(x, y) 6 
M 2 : x > 0,y > 0}, except the one corresponding to the unique fixed point, are the ovals of 
some elliptic curves of the form 

{V a (x, y) = h} = {(x + l){y + l)(x + y + a) - hxy = 0}. (2) 

The action of the map on each of the above curves can be described in terms of the group 
law of them: in fact, the Lyness' one is a particular case of the well known family of QRT 
maps, In particular all possible periods of the recurrences generated by ([I]) are known, 
and for any a ^ {0, 1} infinitely many different prime periods appear (the cases a = and 
1 are globally periodic with periods 6 and 5 respectively). 

Recently, there has been some progress concerning the study of the non-autonomous 
periodic Lyness' equations 

a n + Un+l ,„x 

U n +2 = , (3) 

U n 

when {a n } n is a /c-periodic sequence taking positive values, and the initial conditions u%, U2 
are, as well, positive (see [6j [TJ [H [181 EL and also [TTJ, [22]). These works focus on some 
qualitative aspects of the dynamics like persistence, stability, as well as some integrability 
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issues. In this paper we will focus on the characterization of periodic solutions of the 
2-periodic case given by 

f a for n = 2!+l. 
a n = { 4 
{ b for n = 21, 

where I £ N and a > 0, b > 0. 

The solutions of equation ^ can be studied through the dynamics given by the com- 
position map 

*W*> V) := (F b o F a )( X , y) = (^±1, a ^±A) , (5 ) 

\ x xy J 

since the relation between the terms of recurrence (^-Q and the iterates of the composition 
map is given by 

(u 2n +l,U2n+2) = F bta (u2n-l,U2n), and (u 2n +2,U2n+3) = F a ^{u 2m U 2n +l), (6) 

where {ui,u 2 ) £ Q + and n > 1. 

The map F b ^ a is integrable in the sense that it has a first integral, [18j . Here we use the 
following one: 

(bx + a){ay + b)(ax + by + ab) 

Vb,a{x,y) = • (7) 

xy 

This means that each map F bA preserves a foliation of the plane given by cubic curves which 
are, generically, elliptic. The composition maps F^ a associated to the 2-periodic Lyness' 
equations are also particular cases of QRT maps. It is interesting, however, to notice that 
the 2-periodic case is one of the few rationally integrable or even meromorphically integrable 
fc-periodic Lyness' equations (see [8]). The maps F^ a are also a generalization of the 
map considered in |14j . 

It is known that a map F^ a has a unique fixed point (x c , y c ) £ Q + given by the solution 
of the system 

[ y A = b + x. 

which corresponds to the unique global minimum of V& in Q + . Furthermore, setting 
h c := {Vi, : a(xc, J/c)} ; the level sets Cf '■= {{Vb A = h} n Q + for h > h c } are closed curves 
and the dynamics of Fj, a restricted to these sets is conjugate to a rotation on the unit circle 
with associated rotation number 0b, a (h), [§]• In Section [2] we give an alternative proof of 
this fact (see Corollary [6]). In this paper we study some properties of the rotation number in 
order to characterize the possible periods that can appear in the family of recurrences Q. 

The main results of the paper are the following. 
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Theorem 1. Set 

I(a,b) := (<j(a,b), ^ 

where (c,d) = (min(c, d), max(c, d)), and 

a + bx c 

1 - - 



cr(a, b) : = — arccos , 



+ x c 



1 

2^ 





1 


arccos - 


-2 + 


V2 


Xc Vc 



For any fixed a,b > 0, and any value 9 G I (a, b), there exist at least an oval of the form 
such that the map restricted to the this oval is conjugate to a rotation, with a rotation 
number 9b, a (h) = 9. 

Notice that F\ \ is the doubling of the well-known globally autonomous Lyness' map F\ 
which is globally 5-periodic and whose rotation number function takes the constant value 
l/fQ an d therefore 1(1, 1) = {2/5}. Notice also, that there are other values of a and b for 



which I (a, b) = {2/5}, as can be seen in Proposition 18 



Theorem 2. Consider the family of maps Fb ja given in (5p fora,b>0. 

(i) If (a, b) ^ (1, 1), then there exists a value po(a, b) G N such that for any p > po(a, b) 
there exist at least a continuum of initial conditions in Q + ( an oval Cf ) giving rise 
to p-periodic orbits of Fb^. Moreover, when o~(a,b) ^ 2/5 the number po(a,b) is 
computable. 

(ii) For each number 9 in (1/3, 1/2) there exists some a > and b > and at least an oval 
C£, such that the action of Fb A restricted to this oval is conjugate to a rotation with 
rotation number 9b A (h) = 9. In particular, for all the irreducible rational numbers 
q/p G (1/3, 1/2), there exist periodic orbits of Fb A of prime period p. 

(Hi) The set of periods arising in the family {F^^, a > 0,6 > 0} restricted to Q + contains 
all prime periods except 2, 3, 4, 6 and 10. 

In fact, as we will see in Section [4j the prime periods 2 and 3 do not appear for any a 
and b in the whole domain of definition of the dynamical system defined by Fb a , but periods 
4, 6, and 10 appear for some a, b > and some initial conditions in Q \ Q + . 

Corollary 3. Consider the 2-periodic Lyness' recurrence (0j-(gp for a > 0, b > and 

positive initial conditions ui and ui- 



'Here we adopt the notation in [I] and [2f5] concerning the determination of the rotation number. Notice 
that in [B^ it is adopted the determination given by 1 — 0b ta (h). 
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(i) If (a, b) ^ (1, 1), then there exists a value p$(a, b) E N such that for any p > po(a, b) 
there exist continua of initial conditions giving rise to 2p-periodic sequences. More- 
over, when o~(a,b) ^ 2/5 the number po (a, b) is computable. 

(ii) The set of prime periods arising when (a, b) 6 (0,oo) 2 and positive initial conditions 
are considered contains all the even numbers except 4, 6, 8, 12 and 20. If a^b, then 
it does not appear any odd period, except 1. 

With our tools, the number po(a, b) in the statement (i) of both Theorem [2] and Corollary 
[3] is only computable if o~(a,b) 7^ 2/5. As can be seen from Proposition [18] or Corollary 19 



this happens in an open an dense subset of the parameter space, and therefore po(a, b) is a 
generically computable value. 

Observe that Theorem [2] characterizes all the prime periods that can appear for the 
iterations of the maps F^ a in Q + . In fact it characterizes the set of periods of any planar 
system of first order difference equations associated to any map conjugate to Fb >a , like the 
two families of first order systems of difference equations 

J u n+ iu n = a + v n , ^ \ u n+ iu n = a(l + v n ), 
[ v n+1 v n = b + Un+i, \ v n+l v n = /3(1 + u n+ i). 

for some choices of a, b,a,(3 > 0, and u\ > 0, v\ > 0. The second one is associated to the 



map G, defined in (11) below. 

The paper is structured as follows. In Section [2] we see that on the level sets C^, the 
action of the map F\, a can be seen as a linear action of a birational map on an elliptic 
curve, and thus we reobtain that on these level sets it is conjugate to a rotation (these 
also follows from the fact of being a QRT map, but the notation and elements introduced 
here will be useful in our further analysis of the rotation number function). The main 
part of this section is devoted to proof the elliptic nature of the invariant curves on M 2 , 
C-h '■= {Vb, a = h}, as well as to derive a Weierstrass normal form representation of both 
the invariant curves and the map. This one, will be a key step in our approach to the 
asymptotic behavior of the rotation number function when h tends to infinity. This is done 
in Section [3j which is devoted to study the limit of the rotation function. In particular, we 
obtain that the asymptotic behavior of this function at the energy level corresponding to 
the fixed point and at the infinity do not coincide. It is worth noticing that further tools to 
study the rotation intervals have to be developed to face the problem of finding the set of 
periods when these limits coincide. The main results, as well as some other ones concerning 
periodic orbits are proved in Section [4} 
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2 Ff) >a as a linear action in terms of the group law of a cubic 

2.1 An overview from an algebraic geometric viewpoint 

In this section we will see Fb >a as a linear action in terms of the group law of the cubic. The 
level sets {Vb ja = h} are given by the cubic curves in M 2 : 

Ch = {(bx + a) (ay + b) (ax + by + ab) — hxy = 0}. 

These curves Ch, in homogeneous coordinates [x : y : t] € CP 2 , write as 

Ch = {(bx + at) (ay + bt)(ax + by + abt) — hxyt = 0}. 

Observe that there are three infinite points at infinity which are common to all the above 
curves 

H = [1 : : 0]; V = [0 : 1 : 0]; D = [b : -a : 0]. 

Notice that none of these points is an inflection point. 
The map Fb A extends naturally to CP 2 as 

Fb,a ([x : y : t]) = [ayt + y 2 : at 2 + bxt + yt : xy] , 

and leaves the curves Ch invariant. 

Notice that only a finite number curves Ch are singular, and thus correspond to non- 



elliptic curves Ch, since the discriminant is a polynomial in h (see Section 2.2.3). In partic- 



ular, we will prove that for all the energy levels h > h c , the curves Ch are non-singular. 
Proposition 4. If a > and b > 0, and for all h > h c , the curves Ch are elliptic. 



The above result will be a direct consequence of Proposition [9] (proved in Section 2.2) 



and it implies that in Q + \ {(x c ,y c )} the map Fb A is a birational transformation on an 
elliptic curve, and therefore it can be expressed as a linear action in terms of the group law 
of the curve ([19, Theorem 3]). Indeed, for those cases such that Ch is elliptic, if we consider 
Po = [^o : Ho '■ 1] ^ Ch, and taking the infinite point V as the zero element of Ch, then 
Pq + H can be computed as follows: take the horizontal line passing through Po and H, it 
cuts Ch in a second point (Pq * H) = (x\,yo) (where M * N denotes the point where the 



line MN cuts again the cubic). A computation shows that x\ = (a + yo)/xQ. The vertical 
line passing through Pq * H cuts Ch at a new point Pi = [x± : y± : 1] = Pq + H . Again a 
computation shows that y\ = (a + bxo + yo)/(xoUo)- Hence (xi,y\) = Fb,a(xo,Uo), and we 
get the following result: 

Proposition 5. For each value of h such that Ch is an elliptic curve, Fb a \~ (P) = P + H , 
where + is the addition of the group law ofCh taking the infinite point V as the zero element. 
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So we have the following immediate corollary 

Corollary 6. On each level set the map Fj, >a is conjugate to a rotation. 

Remark 7. Notice that the zero element V of the group law is not an inflection point 
of the cubic. So the relation of collinearity for three points A, B, C does not mean that 
A + B + C = V , but A + B + C = V *V (an inflection point U satisfies U *U = U ). Another 
useful relation is —P = P * (V * V), where V * V = R := [—a/b : : 1]. 

The above results are a consequence of the ones proved in the next section. Now we are 
able to give a first result about admissible periodic orbits of F^ a . 

Proposition 8. For any admissible period p of F^ a , there is only a finite number of level 
curves of the admissible invariant levels sets Ch, which may have this period. 

Proof. As we have noticed, there is only a finite number of energy levels {Vb, a = h}, which 
correspond to non-elliptic curves. For the rest of values (those such that Ch is an elliptic 
curve) we have the following: if there is a periodic orbit with minimal period 2n in the 
curve Ch D Q + with h > h c , since F^ a is conjugate to a rotation, then Ch must be filled 
by with periodic orbits of minimal period 2n, and also, taking into account the natural 
extension F^ a (defined in Proposition [HJ we have 2nH = V or, in other words, if and only 
if nH = —nH. Using Remark [7] we have that —nH = nH * R, where R = (—a/b, 0, 1). So 
there exists 2n-periodic orbits if and only if nH = nH * R, that is, if and only if R is on the 
tangent line to the curve Ch at the point nH. Computing this tangent line and imposing 
this condition, we obtain that the energy level of corresponding to this curve must be a zero 
of a polynomial P2n(h), so it has a finite number of solutions. 

Suppose that Ch is filled by periodic orbits of period p = 2n + 1, then 

(2n + l)H = V ^ nH + H = -nH ^ (nH * H)*V = nH * R, 

which is again an algebraic relation on h giving rise to a finite number of solutions. | 

2.2 Proof of the ellipticity of the curves Ch for h > h c and a Weierstrass 
normal form representation for F& >a 

In this section we prove Proposition |4j In order to simplify the computations we introduce 
a preliminary change of variables and parameters. It is easy to check from ([T]) that the 
integral V\, a admits the following form 

w,)/m=Ci + ae^)(Hj + D/[(i)®]. 
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This means that taking the new variables 



X:=%, Y:=y, a: =«p:=^, (9) 
b a b z a z 



the first integral is now given by 



the map F^ ta becomes 



w{XY)= (X + c,)(Y + m + X + Y) ; 

X Y 



G(A -,V) = (2fl±H, ^ + A ° y + °) ); (1!) 



and the curves Ch are brought into 

V L = {(X + aT)(Y + f3T)(T + X + Y) - LXYT = 0}, 

where L = h/(a 2 b 2 ). Observe that all the curves T>l have three common infinite points 
given by H = [1 : : 0]; V = [0 : 1 : 0]; and D = [1 : — 1 : 0], and again, none of these points 
is an inflection point. The map G extends to CP 2 as 

G ([X : Y : T]) = [a(YT + Y 2 ) : f3(XT + aYT + aT 2 ) : XY] 

The level L c = h c / (a 2 b 2 ) corresponds to the level of the global minimum of the invariant 
W in Q + . Proposition |4] is, then, a direct consequence of the following result, that will be 
proved in Subsection |2.2.3| 

Proposition 9. If a > and f3 > 0, for all L > L c = h c /(a 2 b 2 ), the curves T>l are elliptic. 

Reasoning as in the case of Fb A we obtain the following result. 

Proposition 10. For each L such that T>l is an elliptic curve, G\ v (P) = P + H , where 
+ is the addition of the group law of T>l taking the infinite point V as the zero elemen^ 

As a consequence of the Propositions [9] and 10 we obtain that the action of G on each 
level T>l n Q + , with L > L c is conjugate to a rotation. 

In order to prove Proposition [9] as well as to significatively simplify some of the proofs 
of the results concerning the behavior of 9b,a(h) done in Section [3j we first present a new 
conjugation between the linear action of G on T>l (and therefore for F^ a on Ch) and a linear 
action on a standard Weierstrass normal form of the invariant curves. We prove: 

' Notice that, again, the zero element V of the group law is not an inflection point of the cubic (see 
Remark [7]). 
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Proposition 11. (i) For any fixed L, the curves T>l have an associated Weierstass form 
given by 

£ L = {[x:y: t], y 2 t = 4x 3 - g 2 xt 2 - g 3 t 3 }, (12) 



being 



92 192 



l - (l 8 + ][>(«, /3)z/) and g, = ^ (^-L 12 + £ ft(a, 0)L*j , 



where 



p 7 (a,b)= -4(a + /3 + l), 

p 6 (a,6)= 2 (3(a - /3) 2 + 2(a + /3) + 3) , 

p 5 (M) = -4(a + /3- 1) (a 2 -4/3a + /3 2 - l) , 

p 4 (o,6)= (a + /3-l) 4 . 



and 



gn(a,6)= 6 (a + + 1), 

gi0 (a, 5) = 3 (-5a 2 + 2a/? - 5/3 2 - 6a - 6/3 - 5) 

©(a, fe) = 4 (5a 3 - 12a 2 /3 - 12a/3 2 + 5/3 3 + 3a 2 - 3a/3 + 3/3 2 + 3a + 3/3 + 5) 

g&(a, b) = 3 (-5a 4 + 16a 3 /3 - 30a 2 /3 2 + 16a/3 3 - 5/3 4 + 4a 3 

-12a 2 /3 - 12a/3 2 + 4/3 3 + 2a 2 - 8a/3 + 2/3 2 + 4a + 4/3 - 5) 

g 7 (a, 5) = 6 (a 2 - 4a/3 + /3 2 - l) (a + /3 - l) 3 

q 6 (a,b)= -(a + /3-l) 6 

(m^ For any L > L c , the curve £l is an elliptic curve. 
(Hi) For any L > L c , G\d l is conjugate to the linear action 

G 1el :P^P*H, (13) 

where 



H 



- 1 - (L 2 - 2 (a + /3 + 1) L + (a + /3 - l) 2 ) L 2 : -^a/3L 4 : 1 
48 V / 8 



(14) 



and-k denotes the sum operation of the cubic E^, taking the infinity point V =[0:1:0] 
as the zero of the group law. 

In the next subsections we give the proof of the above result. 
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2.2.1 Computing the Weierstrass normal form of T>l 

In this section we explicit the transformations which pass T>l into their Weierstrass normal 
form. Although it is a standard calculation, we prefer to include it because in the proof of 



Proposition 1 1 , we need the explicit expressions of some of the intermediate transformations 
below. 

First transformation: Following |25| Section 1.3], we take the following reference projective 
system: Let t = be the tangent line to T>l at V. This tangent line intersects T>l at 
R = [—a : : 1], so we will take the reference x = to be the tangent line at R. Finally 
y = will be the Y axis. The new coordinates x, y, t are, then, obtained via 

x = P(l- a)X + aLY + a(3(l - a)T; y = X;t = X + aT, (15) 

and the equation of the cubic becomes 

xy 2 + y(a\xt + a 2 t 2 ) = a^x 2 t + a^xt 2 + a^t 3 , (16) 

where 



a x = {a-p- 1-L)/L, a 2 = P ((1 - a) 2 + /3(1 - a) - L) /L, a 3 = -1/L 



2 



a 4 = (2/3(1 - a) - L(l + p)) /L 2 , a 5 = P(p(l - a) - L)(L - 1 + a)/L 2 



(17) 

Second transformation: We multiply both members of ( |16| ) by x, and we take the new 
variables (denoted again X, Y, T): 

X = x;Y = xy/t;T = t, (18) 

obtaining the new equation 

T 2 Y 2 + ( ai X + a 2 T)YT 2 = a 3 X 3 T + a 4 X 2 T 2 + a 5 XT 3 . (19) 



Third transformation: Removing the factor T in (19) and taking the new variables x,y,t 

x = X;y = Y+{a 1 X + a 2 T)/2;t = T, (20) 
we obtain the new equation 

ty 2 = Ax 3 + [ix 2 1 + vxt 2 + 7t 3 , 

where 

A = 03, u. = 04 + a\ /4, v = 05 + 0102/2, 7 = a|/4. (21) 
The last two transformations Let X, Y and T be the new variables given by 

X = x/(4X), Y = y/(4X 2 ), T = t (22) 
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In these new variables the curve writes as 



TY 1 



Finally, with the last change 



A 2 4A 3 4 2 A 4 



rpd 



1 12X\ 2 + nt 



12 A 2 

the curve is written in the Weierstrass' normal form 



y = Y, t 



(23) 



y 2 t 



4x s - g 2 xt 



93t 3 , 



where 



92 



12A 4 4A 3 



93 



+ 



7 



216A 6 48A 5 16A 4 ' 



(24) 



2.2.2 Proof of Proposition 11 



Proof of Proposition 11 (i) The expression of the equation ( 12 ) comes from the transforma- 



tions leading T>l to its Weierstrass form given in Section 2.2.1. In particular the expressions 



of the coefficients of g 2 and g% in terms of a, /3 and L come from formulae (17), (21) and 



<|24|). 

(ii) From Proposition [9j the curves T>l, with L > L c are elliptic. This means that for 
L > L c the curves £l will be also elliptic. 



Statement (iii) is a direct consequence of Proposition 10 But it is still necessary to 
prove that the zero point V is transformed into V = [0 : 1 : 0]. To do this we should obtain 
the image of V under the same transformations as above, but the image of V under the 
first transformation, say V\ := [1 : : 0], has no image for the second transformation. So 
we use a continuity argument: Set V\{x) = [x : y(x) : 1], where y(x) is choose in sort that 
Vi(x) is on the cubic (16), that is, y(x) can be either y+{x) or y~(x), where 



y±{x) 



-a\x — a>2 ± 4 CL3X 3 + 4 (04 + ai 2 )x 2 + 2 {a\a 2 + 4 as) x + a 2 2 

2x 



The leading term of y±(x) when x — > +00 is ±^/a3 x (observe that these two numbers are 
imaginary, because 03 = — 1/L 2 ). So we have two complex points of the cubic [x : y±(x) : 1] 
where y± = ±^/a3 x + 0(1): 

Vi(x) = [x : ±^3~x~ + O(l) : 1] = [1 : + 0(l/x) : 1/x] 

which tends to V\ = [1 : : 0] when x tends to +00. 

But if we take the images of these two complex points, in their first form, by the 



transformation (18), we obtain 

V 2 (x) = : ±v^3 + 0(l/y/x) : 1/v 7 ^], 
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which tends to [0 : ±y/ai : 0] = [0 : 1 : 0]. So this point is V2, the image of V\ obtained by 
continuity. Now the successive transformations can be done without any problem and we 
obtain V = [0 : 1 : 0]. 



The point H is the image of H = [1 : : 0] under the successive transformations (15) 



(18), (20), (22) and (23). | 



2.2.3 Proof of the ellipticity of the curves T>i for L > L c 

In this section we prove Proposition [9j but first we notice the following facts: 

Fact 1. The curves T>l have no singular points at the infinity of CP 2 . Indeed, setting 
D := (X + oT){Y + f3T)(T + X + Y) - LXYT, we have that 

|^(X,y,0) = Y(Y + 2X), ^(X, Y, 0) = X(2Y + X), 

so an straightforward computation shows that on the line T = the only singular point 
must satisfy X = and Y = 0, thus not in CP 2 . 

Fact 2. An straightforward computation shows that the set of affine singular points of T>l 



coincides with the set of singular points of the invariant function W given in (10). On the 
other hand, a computation shows that 

W' x = (Y + (3)(X 2 -a{l + Y))/(X 2 Y) and 



(25) 

— ( \ 4- a II)'-- si 1 -!- A \ \ ' I \ y-l 

y 



W' ={X + a)(Y 2 - p(l + X))/{XY 2 



So for L > L c , the other possible singular point of a curve T>l is one of the point intersection 
of the two parabolas Vi := {X 2 = a(l + Y)} and V 2 := {Y 2 = (3(1 + X)}. 
Fact 3. The equilibrium point P = (p, q) which is in the positive quadrant, is one of the 
(at most) four real points on V\ CiT , 2, and it is on T>l and singular in this curve only for 
L = L c (in this case, it is a real isolated point). 

Now we introduce the some new parameters (p, q) which correspond to the coordinates 
of the fixed point in Q + , thus given by 

f *> = «(! + .), (26) 

\ 9 2 = «l+p). 

With these new parameters it is easy to obtain an explicit form for L c : 

_ (1+p + q) 3 
Lc ~ (l+ P )(l+9)- {27 > 

In fact this explicit expression allows to obtain an straightforward proof that L c > 1 
since 

L c =( P + ,+ l)(l + j ^ T ) (l + ^i)>l- 
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Furthermore, after some computations (skipped here) it is possible to prove that 



lim >,g) = (0,0). 
(a,/3)->-(0,0) 

This implies that the above one is the sharpest lower bound for L c . 

Fact 4. Any singular point ofT>L with L > L c must be real. If a singular point is not real, 
the conjugate point is also singular and different, and so the (real) line which support these 
two points is in the curve, which splits in a right line and a conic. It is easy to see that, 
after the identification of the coefficients the equation of the cubic must be written as 

D = (X + Y + (1 - L)T) (xY + j3XT + aYT + y^jT 2 ^ = °- 

The point [0 : —1 : 1] is on the cubic, hence — L (—a + a/3/(l — L)) = 0, which is impossible 
because the above number is positive, since we have just seen that L > L c > 1. 

Taking into account all the above information, to prove Proposition [9j we only need to 
show that if T>l is singular then L < L c . 

Lemma 12. // the curve T>l is singular, then L < L c . 



Proof. Using the Weierstrass form given by (12) and the expressions of 52 and 53 given by 



from formulae (17), (21) and (24). We obtain that the discriminant of T>l is given by 



A := g 3 2 - 27<? 3 2 = ^a 2 f3 2 L 15 A 1 (L;a, (3), 

where Ai(L) = L 4 + Yli=o <>i( a i P) L l , and being 

5 3 (a,f3)= af3 - 4a - 4/3 - 4, 

S 2 ( a , (3) = -3a 2 (3 - 3af3 2 + 6a 2 - 21a/3 + 6/3 2 + 4a + 4/3 + 6, 

ft (a, /3) = 3a 3 /3 - 21a 2 /3 2 + 3a/3 3 - 4a 3 + 18a 2 (3 + 18a/3 2 - 4/3 3 + 4a 2 

-25a(3 + 4/3 2 + 4a + 4/3 - 4, 
5 (a,/3)= -(a-l)(/3-l)(a + /3-l) 3 . 

Obviously L = gives a singular curve (a union of three lines). To see that all the real 
roots of Ai are lower than L c , we take into account the explicit expression of L c given by 



(27), we use the change 



P a Q 

1 + q l+p 

and also the translation given by L = L + L c , obtaining that Ai(L; a, /3) = if and only if 
the numerator of Ai(L + L c ;p 2 /(1 + q),q 2 /(l + p)) vanishes, which happens if and only if 
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At(L):=L £? =0 di(p, g )L' 



0, where 



d2(p,q) 
di(p,q) 



(l+p) 3 (l + g) 3 , 

(1 + p) 2 (l + q) 2 (p 2 q 2 + I2p 2 q + 12pg 2 + 8p 2 + 20pq + 8g 2 + 8p + 8q) , 
(1 + p) 2 (l + g) 2 (9pV + 9p 3 q 4 + 60pV + 93pV + 60p 2 g 4 + 64p 4 g + 228p 3 g 2 
+228p 2 q 3 + 64pg 4 + 16p 4 + 176p 3 g + 312p 2 g 2 + I76pq 3 + 16g 4 + 32p 3 + 160p 2 g 
+ I60pq 2 + 32g 3 + 16p 2 + 48pq + 16g 2 ) , 
cfo(p, q) = pq(p 2 +pq + q 2 +p + q) (3pq + 4p + 4g + 4) 3 . 

Since Ai is a polynomial with positive coefficients, its real roots are negative or zero, and 
thus, the real roots of Ai are less or equal than L c . | 



3 Behavior of 6^ a (h). 

To prove Theorem [TJ as well as to study the periods appearing in the family F^ a we need 
to analyze the rotation number function 0b a (h), with special emphasis on its asymptotic 
behavior when h — > hj, and when h — > +oo. Regarding this last question, in |6j, there 
were pointed out strong numerical evidences that lim 9b, a (h) = 2/5. Notice that these 
evidences agree with the fact, proved in p], that for the autonomous case the rotation 
number associated to the invariants curves (2), say 9 a (h), satisfies lim 9 a (h) = 1/5 so it 

h— >+oo 

is independent of the value of the parameter a. In this section we prove 
Proposition 13. The following statements hold 

(i) The rotation number map 9ba(h) is analytic on (h c ,+oo). 



(ii) lim 0b, a (h) 



(Hi) The rotation number map 0b a (li) is continuous in [/i c ,-t-oo). Furthermore, setting 



it holds that 



cr(a, b) 



a(a,b) := lim 6 b ^ a (h), 

h^ht 



1 (\ \ a + bx c 

— arccos - 
2vr I 2 



1 f 1 

2^ arCC ° S 2 



-2 + 



x c (b + x c 
1 



%c Uc 



(28) 

(29) 
(30) 



To prove the statements (i) and (ii) of Proposition 13, it will be helpful to use the 



conjugation of Fb a given by the map G defined in (13), for which its invariant level sets are 
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described by the Weierstrass normal form Sl denned in (12), and then to apply a similar 
machinery as in pp. First, we recall that given the elliptic curve £l, for L > L c , there exist 
two positive numbers lj% and u; 2 depending on a, (3 and L and a lattice in C 

A = {2nu!\ + 2miui2 such that (n, m) E Z 2 } C C, 

such that the Weierstrass p function relative to A 

p{ z ) 



aga\{o} 



1 



1 



>-A) 2 A 2 

gives a parametrization of fx- This is because the map 

0: C/A — > £ L 

_^ f : p'(*) : 1] if z$ A, (31) 

\ [0 : 1 : 0] = V" if z G A, 

is an holomorphic homeomorphism, and therefore 

p'(z) 2 = Ap{zf - g 2 p(z)-g 3 (32) 

(see [20, Lemma 5.17] and [24, Proposition 3.6]). Another interesting property of the 
parametrization of £l given by <f> is that it is injective for the real interval (0, ujx) onto the 
real unbounded branch of £l whose points have negative y -coordinates, and that <j>(u)\) = 
[ei, 0, 1] or, in other words: 

ei = p(ui) (33) 
(see again [201 El])- Finally, observe that (f)(0) = <f>(2ui) = V, and lim p(u) = +oo. Taking 

u— >0 

all the preceding considerations into account, and by direct integration of the differential 



equation (32) on [0, u), we have that in real variables 

-+oo 



U 



ds 



92S ~ 53 



(34) 



Proof of Proposition |i3| (i). Observe that equation (12) can be written as 



where 



ty 2 = 4(x — te\)(x — te2)(x — tes), 



ei + e 2 + e 3 = 0, e x e 2 + e 2 e 3 + e 3 ei = -g 2 /4, eie 2 e 3 = g 3 /4, 



and e 3 < e 2 < e\. Then, equation (34) writes: 



u 



ds 



P (u) y/4(s - ei)(s - e 2 )(s - e 3 ) 



(35) 
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and from the above relation and equation (33) we get: 

" + 0O 



ds 



■l V^( S - e l)( S - e 2)( s - e 3) 



From Proposition 
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the action of G on each level £l, with L > L c is conjugate to a 
rotation of angle 2tt@(L), where is the rotation number associated to G\g L - On the other 
hand from the expression of H given in (14) we know that its y-coordinate is negative, so 
from the above considerations on the parametrization <p we know that there exists u £ (0, ui) 
such that it has the following representation 



/2vri 





exp (2vrie(L)) <G [0, 2uj{\/K = S 1 , 



so u = 2ui@(L), and since H is the image of V under the transformation G\g L , it is given 
by 4>(2ojiQ(L)), and therefore 

p(2wi0(L)) = X{L), 



where X(L) is the abscissa of point H given in (14). Applying (35) we obtain that 

-+oo 



2wi6(£) 



ds 



hence 



2G(L) 



x(L) ^4(5 - ei)(s - e 2 )(s - e 3 ) 

+°° ds 
x(L) V(s ~ ei)(s - e 2 )(s - e 3 ) 
ds 



(36) 



ei)(s - e 2 )(s - e 3 ) 

By using the change of variables s = e\ + 1/r 2 and r-^/ei — e 3 = n, it can be proved that 



equation (36) becomes 



20(L) 



' e l~ e 3 



dn 



V(l + n 2 )(l + en 2 ) 
du 



(37) 



o ^(l + n 2 )(l + eu 2 ) 

where ^ = A(L) — ex, and e = (ei — e 2 )/(ei — e 3 ). 

The analyticity of 0(L) in (L c , +oo) can be derived now from expression (37), by taking 
into account that when L > L c , e\, e 2 , e 3 are different, and that from the geometry of the 



problem we have e 3 < e 2 < e\ < X(L), so all the elements in the left hand side of (37) are 
holomorphic in a neighborhood of (L c , +oo). Hence @(L) is analytic on (L c , +oo). | 



Observe that lim 9b, a {h) 



lim O(L), hence our main objective now is to prove that 

L— ¥+oo 



lim Q(L) = 2/5. To do this we will need the following auxiliary result 
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Lemma 14. ([I, Lemma 4]) Let A, £,7 be positive numbers. For any map (p(e) such that 



lim (j)(e) = 0, and A + 4>(e) > 0, set 

e^O 



iV(e,A, 7 ) 



A+-Ke) 



du 



z, and D{e) 



du 



+00 

7(l + n 2 )(l + en 2 )' 



y/(l + v?)(l + ev?) 

T/ien D(e) ~ (1/2) ln(l/e), and if j < 1/2 we have N(e, A, 7) ~ 7 ln(l/e), where ~ denotes 
the equivalence with the leading term of the asymptotic development at zero. 



Proof of Proposition 
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(ii). Following the steps in pQ, we will compute lim @(L) by 

L— H-oo 



using the asymptotic developments of the elements involved in equation (37). 

Solving the cubic equation 4 s 3 — 52 s — #3 = by leaving the solutions in terms of 52 and 



#3; using the expressions of these coefficients in terms of a, (3 and L given by Proposition 11 



and after some computations done with a symbolic algebra system (Maple 12), we obtain 
that 

ei ~ L 4 /48, e 2 ~ L 4 /48 and e 3 L A /2A, 

when L tends to infinity, and 

ei — e2 ~ af3 L 3 ^ 2 and e\ — e% ~ L 4 /16. 

Setting e = (ei — e2)/(ei — 63), we obtain e ~ 16a/3/L 5//2 (when L tends to infinity), so 

2 2 3 / 5 (a/3) 2 / 5 



L 



r2/5 



(38) 



(when £ tends to zero). From the expression of H given in (14) we have 

v = X(L) - ei ~ a/3L 2 /4, 



and therefore, using (38), we get 



ei - e 3 



L 



A£~ 2 I\ 



where A = 2 3 / 5 (a/3) 1 /( 10 ). So from equation (37), and Lemma 14 we finally obtain 



29(L) 



N(e,A,2/5) |ln(l/e)_4 



D{e) 



iln(l/ £ ) 



and therefore lim @(L) = 2/5. | 
It is well known (see [2] for instance, and also [26]) that under some regularity con- 
ditions the rotation number function of a differentiable planar map F near an elliptic fixed 
point, can be extended continuously to the fixed point. In this sense, the proof of statement 
(ii) follows by applying the following version of [21 Proposition 8]. 
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Lemma 15. Let U be an open set o/R 2 and a C 1 map F : U — )■ U. Suppose that F has 
a unique elliptic fixed point P in U (that is DF(P) has two distinct non real eigenvalues 
with modulus 1). Suppose that F has a first integral V in U, attaining a strict minimum 
hp := V(P) at the point P, and such that 

(i) For h > hp, the invariant level sets {V = h} (except the level set {V{P) = hp}) are 
closed curves, surrounding P. 



(ii) For h > hp, the sets Int ({V = h}) are starlike with respect to P, and each half line 
from P cuts {V = h} n U exactly in one point. 

(Hi) On each curve {V = h} for h > h p , the map Fuy—^x is conjugate to a rotation of 
angle 2ir9{h). 



Then, 



9{P) := lim 9(h) = ±- arccosf ^g^ ) . (39) 
h^hp 2tt V 2 ' 



Proof of Proposition \13\ (Hi). As a direct consequence of Lemma 15 , the continuous extension 
of 9^ a (h) to is 

o-(a,b) = -'-arccos ( -Trace(DF b Jx c , y c )) 



Then, the expressions (29-30) are then easily derived from ([8]) after some easy manipula- 
tions. ■ 



4 Periodic solutions 

4.1 The rotation interval I(a,b). Proof of Theorem [l] 

In this Section we prove Theorem [TJ and we give a more precise characterization of the 
rotation number interval I(a,b). 

Proof of Theorem [IJ In proposition 13, it is proved that 9b, a (h) is an analytic function on 



(h c , +oo) and continuous in [h c , +oo]; that lim 9f )a (h) = a(a, b); and lim 9t, a (h) = 2/5. 

h-*h£ ' h-^+oo 

It is clear then that for any 9 in the interval I(a, b) defined by 

I(a, b) := ( lim 9 b a (h), lim 9 b a (h)\ = (cr(a,b) 
\h->h+ h^+oo I \ 



' r / ' 



where (c,d) = (min(c, d), max(c, d)), there is a value of h G [h c , +oo) such that the action 
of F b;a restricted to Ct is a rotation with rotation number 9, since 

I(a, b) C Image {9 b Mh c , +oo))) . 
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To have a better description of the interval I (a, b), it would be interesting to know the 
relative position of a(a,b) with respect to 2/5 for given values of a and b. This is issue is 
considered in Proposition 18 and Corollary 19 (see also Figure 1), but before proving them 
we need a preliminary result. 

Set J(a, b) := Image (9b, a ((h c , +oo))) C [0, ^]. As as stated above, it holds that I(a, b) C 
J(a,b). Observe that both intervals are the same when the function 6b, a {h) is monotonic, 
but notice that in [6] it is proved that there are values of a and b for which the map is 
non-monotonic. With respect the interval J(a,b), we prove: 

Lemma 16. The following statements hold. 

(i) For any a,b > 0, a (a, b) > —. 

(ii) lim a(b 2 ,b) = -. Therefore, when b is sufficiently large o~(b 2 ,b) > 2/5, and then 
b-^+oo 2 

I(b 2 ,b) = {2/5,a(b 2 ,b)). 

(Hi) lim a(b 2 ,b) = -. Therefore, when b is small enough a(b 2 ,b) < 2/5, and then 
&->o+ 3 

I(b 2 ,b) = {a(b 2 ,b),2/5). 
As a first direct consequence of the above result, we have 
Proposition 17. It holds that 



a>0,fe>0 



Furthermore 



7 6>0 



Proof of Lemma 16 (i) Using that 

<j(a,b) = ^arccos QTrace(DF 6j0 (a; c , y c ))Y 

we have that a(a,b) > 1/3 if and only if Trace (-DF&^Xc, y c ))/2 < cos(2-7r/3) = —1/2 or, 
equivalently using equation ( 29 ) , 

_ 1 _ a + bx c _ 1 

x c (b + x c ) 

which is true, since a, b and x c are positive. 
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(ii) From equation (ph we have that y c > \fb. Suppose that a = b 2 , then x c = yja + y c > 



l~a = b. So lim x c = lim y c = +00 and then, using equation (30) 

6— >+oo b— >+oo 



lim 6(b , b) = lim — arccos 
6— >+oo b— >+oo 2ir 





1 


G 


-2 + 




Xc yc 



1 . . 1 

— arccos(-l) = -. 



(iii) When a = 6 2 , we consider the polynomials R(x, y) := x 2 — y—b 2 and S(x, y) := y 2 —x — b 
whose common zeros are the fixed points of Fi a . Observe that y c must be a root of the 
resultant 

Resultant(i?, S; x) = -y(-y 3 + 2by + 1), 
so y% = 2by c + 1, thus y c > 1 and y c <y\ = 2by c + 1, hence 

yc< T ± 2b <2Xb<\. 



So if b < 1/4 we have 



1< yl = 1 + 2by c <l + 4b; 



since the right hand side of the above inequality tends to 1 when b tends to 0, we obtain 
that y c tends to 1. Since x c = y 2 + b, we obtain that lim x c = lim y c = 1, hence 



b~^o+ 



lim 6(b 2 ,b) = lim — arccos ( - 
6^0+ 6^0+ 2tt x 







G 


[-2+ 1 1 




x c y c _ 



6^0+ 



1 < u 1 

— arccosl— 1) = -. 

2tt v ; 2 



The next result characterizes the relative position of a(a,b) with respect 2/5. 

Proposition 18. Let (x c ,y c ) be the fixed point of the map given by Set R(x,y) := 
(j)y — bx — a. Then, a(a,b) < 2/5 if and only if R(x c ,y c ) > 0; a(a,b) = 2/5 if and only if 
R(x c , Vc) = 0; and a(a, b) > 2/5 if and only if R(x c , y c ) < 0. 



Proof. Observe that using expression (29), the condition a(a,b) < 2/5 gives 



thus 



a + bx c 
x c (b + x c ) 

a + bx c 



.'4vr\ 1 + ^/5 

- cos, tJ = — 4—' 



< 



x c (b + x c ) 

where <f> = (1 + \/5)/2. By using (p| this last inequality turns to 



<p (a + bx c ) <Xc{x c + b) = x 2 + bxc = a + y c + bx c . 
Using again <f) — 1 = l/<f>, we obtain that 0y c — 6x c — a > 0. 
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In order to obtain a characterization of the condition a(a,b) = 2/5 in terms of the 
parameters a and b instead of x c and y c , we isolate a and b from the equations ^ and plug 
them in the equation R(x c , y c ) = 0, obtaining 

y c {4> - x c y c + 1) = o. 

Using that <fi + 1 = (j) 2 , we obtain that R(x c , y c ) = if and only if y c = 4> 2 /x c , thus taking 
the parameter t := x c we have that 



t- 3 



and b 



t t 2 ' 

A simple computation shows that the above formulae describe a curve with a single 

2 4 

branch in the parameter's space V := {(a, 6), a, 6 > 0} when 03 < t < <pa . In summary, we 
have proved the following corollary: 



Corollary 19. The curve a(a, b) = 2/5 for a, b > is given by 
T := {a(a, b) = 2/5, a, 6 > 0} = j (a, b) -- 



t^-6 2 6*-t 3 



,t€ (05,<^)!> CP. 



t ' t 2 

It is easy to check that this curve cuts the axes at the points (0, <r*) and (<r*,0) where 
a* = 05 ~ 2.2300, and it splits the parameter space into two connected components. By 
using Lemma [16] (ii) and (iii), it is easy to obtain that one of them corresponds to the set 
{<r(a, b) < 2/5} and the other to {cr(a, b) > 2/5}, and that they have the geometry depicted 
in Figure 1. 



o(a,b)>2/5 




(a,b)=2/5 



Figure 1: The sets {a(a,b) < 2/5} and {cr(a, b) > 2/5}, and the curve a(a,b) = 2/5. 
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4.2 Forbidden periods and periods not appearing in Q + 

Proposition 20. The map Ff,^ a do not have periodic orbits of prime periods 2 and 3. 

Proof. Observe that the orbits of the dynamical system defined by the iterations of F^a are 
also given by the planar first order system of difference equations given by 

J u n+1 u n = a + v n , 
[ v n+ iv n = b + u n+1 . 

If the above system has a periodic solution of prime period 2, then we have that u n +\u n 
is constant, so it also must be a + v n , and therefore v n and u n are constant. 
If there is a periodic solution of prime period 3 of the above system we have 

U n +2Un+lU n = C = U n+2 (a + V n ), (40) 

where c is a constant; and also we have 

Vn+2V n +lV n = d = (b + U n+2 )v n , 

where d is a constant. Hence 

d h 

U n +2 = b. 



Using equation (40) we have c = (a + v n ) yy — bj , so 

bv\ + (ab + c — d)v n — ad = 0, 
and therefore v n is constant. | 

Proposition 21. The map F{, a do not have periodic orbits of prime periods 4, 6 and 10 in 
Q + . 

Observe, however, that there can exist periodic orbits of minimal period 4, 6 and 10 in 
IR 2 \ Q + for some values of a and b. For instance, if a = 1, b = 2, then the curves with 
h = 3; h = 4 + \/l3; and /i ~ 9.8309187775 gives rise to continua of initial conditions with 
minimal period 4, 6 and 10 respectively. These curves are located in M 2 \ Q + . Of course, 
the period 10 appears when considering the map but not as a prime period. 
Proof of Proposition 21 We consider again the coordinates X, Y and T introduced in ([£]), 
and the curves T>l. Observe that, there exist periodic orbits with minimal period 2n in the 
curve T>l if and only if 2nH = V, or in other words, if and only if nH = —nH. Using that 
V * V = R := [—a : : 1], we have 

nH = V <=^- nH = -nH = nH*(V*V) nH = nH * R. 
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The last of the above condition means that there exist a 2n periodic orbit if and only if R 
belongs to the tangent line to T>l in nH. In the following we will see that for periods 4, 6 
and 10 this condition implies L < L c , thus proving that these periods do not appear in Q + . 

Period 4. Since 2H = (H * H) * V = [0 : —1 : 1] =: Q, there exist a 4— periodic orbit if and 
only if R belongs to the tangent line to T>h in Q. An straightforward computation allows 
us to find the explicit expression of the tangent line to 2H: Y = m(L; a, (3) X + n(L; a, /3), 
where 

m(L;a,/3) = — — — r— and n(L;a,/3) = —1. 

a(p — 1) 

This tangent line passes through R if and only if 

— a m(L; a, (3) + n(L; a, (3) = 0, (41) 

or equivalently, after some manipulations, if and only if L + (/3 — 1) (a — 1) = 0. That is, 
when 

L = L 4 := - (P - 1) (a - 1) . 

To see that L4 < L c , we consider again the parameters p and q given by the coordinates 
of the fixed point in Q + , given by (26), that is, we use the change 

l+q l+p 



obtaining that 



(-q 2 +p+l) (p 2 -q-l) (l+p + q) 3 
4 " (l+p) (l + «) < (l+p)(l + q) 



The last inequality can be obtained by subtraction of the numerators, for instance. 
Period 6. A computation shows that 

3H = [-L - a (0 - 1) : LP + p (j3 - 1) (a - 1) : - 1] 

Proceeding as in the period 4 case, we compute the tangent line to T>l at 3H, and we im- 



pose condition (41 ), obtaining that there exists a continuum of periodic orbits characterized 
by the energy level L if and only if either L = < L c or if L is a root of 

P 6 (L) := L 2 + (3aP - la - 2/3 + 1) L + (p - 1) (a - 1) (2q/3 - a - p + 1) . 



Using the change (42) and the translation L = L + L c we obtain that P&{L) = if and only 
if L is not a root of 

P 6 (L) := (1 + q) (1 + p) L 2 + (3p 2 q 2 + 6p 2 q + 6pq 2 + 4p 2 + Upq + Aq 2 + 7p + lq + 3)L 
+ (1 + qf (1 + pf (2pq + 3p + 3q + 3). 
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Observe that all the coefficients of Pq(L) are positive when p and q are positive, so all its 
real roots are negative. Therefore, all the real roots of Pq are less than L c . 
Period 10. A computation done with the aid of a computer algebra system gives 5H := [A : 
B : C], where 

A:= - [{a/3- 1) L + a (3 (/3 - 1) (a - 1)] \l 2 + a (/3 - 1) (/3 - 2) L - a ((3 - l) 3 (a - 1) 
B := -(3 [L + ((3 - 1) (a - 1)] [L + a((3 — 1)] [L 2 + (-2a + 1 - 2/3 + 3a/3) L+ 

+ 03 - 1) (a - 1) (2a/3 - a - /3 + 1)] , 
C := [PL + ((3-1) (J3(a - 1) + 1)] [L + (J3- 1) (a - 1)] [(a/3 - 1) L + a/3 (/3 - 1) (a - T 



We compute the tangent line to V>l at 5-ff and we impose condition (41), obtaining that 
there exists a continuum of periodic orbits characterized by the energy level L if and only 
if 

P W) a(L; a, (3) P 10tb (L; a, (3) P WtC (L; a, (3) = 0, 

where 

P w , a (L; a, /3) = (a /3 - 1) L + a /3 (/3 - 1) (-1 + a) , Pi , b (£, a, /3) = {L + a(/3 - l)) 2 ; 

and Pio, c (^; a, /3) := Ei=o^( a ' Z 3 ) ^> bein g 
p 5 (a,/3)= a^ + 1, 

p 4 (a, /3) = 5a 2 /3 2 - 5a 2 /3 - 5a/3 2 + 12a/3 - 4a - 4/3 + 1, 

p 3 (a, /?) = 10a 3 /3 3 - 20a 3 /3 2 - 20a 2 /3 3 + 10a 3 /3 + 56a 2 /3 2 + 10a/3 3 - 41a 2 /3 - 41a/3 2 

+6a 2 + 35a/3 + 6(3 2 - 6a - 6/3 + 1, 
p 2 ( a , /3) = (/J - 1) (« - 1) (l0a 3 /3 3 - 20a 3 /3 2 - 20a 2 /3 3 + 10a 3 /3 + 55a 2 /3 2 + 10a/3 3 

-37a 2 (3 - 37a/3 2 + 4a 2 + 30a/3 + 4/3 2 - 5a - 5/3 + l) , 
Pl (a, /3) = (/3 (a - 1) - l) 4 (5a 3 /3 3 - 10a 3 /3 2 - 10a 2 /3 3 + 5a 3 /? + 25a 2 ^ 2 + 5a/3 3 

-15a 2 /3 - 15a/3 2 + a 2 + 12a/3 + /3 2 - 2a - 2/3 + l) , 
p (a,/3)= a/3 (/3 - l) 5 (a - l) 5 . 

It is easy to check, using the expression of 5H given above, that when a/3 -1/0 the 
value of L corresponding to the root of Pio,a give 5H = V ', thus characterizing orbits of 
prime period 5, and so 10 is not a prime period for the orbits in this level set. If a/3 — 1 = 0, 
then Pio,a does not depends on L, and it is zero only if a = (3 = 1, which corresponds to 
the well-known globally 10 periodic case a = b = 1. 



Observe that, by using (42), the root of Pio,b is 

p 2 (l +p - q 2 



L = a(l-f3)- 



which is obviously less than L c . 
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To prove that the real zeros of P\q )C are lower than L c , we proceed as in the period 6 



case, by using the change (42) and the translation L = L + L c , obtaining that Pio :C (L) = 
if and only if the degree 5 polynomial P\q^ c {L) obtained when considering the equation 
Pio,c(L + L c ) = with parameters p and q, vanishes. Once again it is easy to check that 
this polynomial has positive coefficients when p and q are positive, so its real roots are all 
negative, and therefore the roots of Pio,c are all lower than L c . g 

4.3 Possible periods. Proof of Theorem [2] 

In order to prove both Theorem [2] and Corollary [3] we need the following auxiliary result. 

Lemma 22. If a ^ b, then the recurrence (^-Q) do not have periodic solutions of odd 
period p ^ 1 . 

Proof. Suppose that for u% > 0,u 2 > the sequence {u n } is a periodic solution with period 
p = 2k + 1. Then, by using the relations in ([6]), we have that 

{ui,U 2 ) = (u 2 k+2,U 2 k+3) = Fa{Fb,a) k (U1,U 2 )- (43) 

But on the other hand, from Q and (43) we have: 

F a (u 1 ,u 2 ) = (u 2 ,u 3 ) = (u 2k+3 ,u 2k+A ) = F b (F ayb ) k (u 2 ,u 3 ) 
= F b (F a , b ) k F a { Ul ,u 2 ) = (F b:a ) k+1 ( Ul ,u 2 ) 
= F b (F a (F b , a ) k ( Ul ,u 2 )^ 
= F b (ui,u 2 ) ■ 

Hence F a (ui,u 2 ) = F b (m, u 2 ), which implies that a = b. | 
Also to prove Theorem [2j we need the following auxiliary result which allow us to 
characterize, constructively, which periods arise when the rotation number takes values in 
a given interval. 

Lemma 23. ([5, Theorem 25 and Corollary 26]) Consider an open interval (c,d); denote 
by Pi = 2,P2 = 3,£>3, . . . ,p n , . . . the set of all the prime numbers, ordered following the usual 
order. Also consider the following natural numbers: 

• Let pm+i be the smallest prime number satisfying that p m +i > max(3/ (d — c), 2), 

• Given any prime number p n , 1 < n < m, let s n be the smallest natural number such 
that p^ n > 4/(d - c). 

. Setp:=p^- l p s r l ---p^-\ 

Then, for any r > p there exists an irreducible fraction q/r such that q/r G (c, d). 
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Proof of Theorem^ (i) If a{a,b) ^ 2/5 (i.e. when (a, b) £ F = {a(a,b) = 2/5, a, b > 0}, 



where T is the curve given in Corollary 19 ), then the value po(a, b) can be computed applying 
Lemma [23] to the interval I(a, b) given in Theorem[TJ as in the proof of statement (hi) below. 

If o~(a,b) = 2/5, then observe that (a, b) = (1,1) is the only point in the parameter 
space V = {(a, b), a, b > 0} such that (a, b) £ {T n {a = 6}}. Hence, as a consequence of 



Lemma 22, if (a, 6) 7^ (1, 1) then J(a, 6) = Image (8b ta ((h c , +00))) is a closed interval with 
nonempty interior (i.e. not the single value 2/5), so there exists some value po(a, b) G N such 
that for all p > po(a,b) there exists irreducible rational numbers q/p S J(a,b). Otherwise 
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6b,a(h.) = 2/5 and all the orbits would be 5-periodic, a contradiction with Lemma 

Observe that the computability of po(a, b) is a generic property in V, since V \ T is an 
open and dense subset of it. 



Statement (ii) is a direct consequence of Proposition 17, which implies that for each 
number in (1/3, 1/2) there exists some a, b > and some initial condition in Q + with this 
associated rotation number for i*& a . In particular, for all the irreducible rational numbers 
q/p € (1/3, 1/2) there are positive values of a and b such that F^ a has continuum of periodic 
orbits of period p. 

we 



(hi) Setting c = 1/3 and d = 1/2, and using the notation introduced in Lemma 23 



have that, m = 7; p\ = 2 (with s\ = 5); P2 = 3 (with S2 = 3); p^ = 5; p$ = 7; ps = 11; 



Pq = 13 and pj = 17 (with Sj = 1 for i = 3, ... ,7). From Lemma 23, for all p G N, such 
that p > po := 2 4 • 3 3 • 5 • 7 • 11 • 13 • 17 = 12 252 240 there exists an irreducible fraction 



q/p 6 (1/3, 1/2). Hence by Proposition 17 there exists some a, b > such that there exists 
a continuum of initial conditions with rotation number 9b a (h) = q/p, thus giving rise to 
p-periodic orbits of F^ a . Now, using a finite algorithm we can determine which irreducible 
fractions q/p with p < po are in (1/2, 1/3), resulting that there appear irreducible fractions 
with all the denominators except 2, 3, 4, 6 and 10. 



Finally, the periods 2, 3, 4, 6 and 10 do not appear as a consequence of Propositions 20 



and [21] | 

Notice that, as a consequence of the arguments in the proof of Theorem [2] (i), we have 
the following result: 

Proposition 24. Set a,b > 0. If (a,b) 7^ (1,1) and <j(a,b) = 2/5, then 9b, a {h) is not 
monotonia on [h c , +00). 

Now the proof of Corollary [3] is straightforward. 

Proof of Corollary [3[ First observe that from Lemma [22| if a 7^ b, then the recurrence 
Q-Q do not have periodic solutions of odd period p 7^ 1. 
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Finally, observe that any ^-periodic orbit of Ff,^ with initial conditions (u\, 112) will give 
rise to a 2p-periodic solution of the recurrence ([3]). Now, both statements (i) and (ii) follow 
as a straightforward application of Theorem [2j | 
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